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Abstract — We present a novel analytical framework for the 
evaluation of important second order statistical parameters, as 
the level crossing rate (LCR) and the average fade duration 
(AFD) of the amplify-and-forward multihop Rayleigh fading 
channel. More specifically, motivated by the fact that this channel 
is a cascaded one, which can be modelled as the product of 
A'^ fading amplitudes, we derive novel analytical expressions 
for the average LCR and AFD of the product of A'^ Rayleigh 
fading envelopes, or of the recently so-called A'^*Rayleigh channel. 
Furthermore, we derive simple and efficient closed-form approxi- 
mations to the aforementioned parameters, using the multivariate 
Laplace approximation theorem. It is shown that our general 
results reduce to the specific dual-hop case, previously published. 
Numerical and computer simulation examples verify the accuracy 
ofHthe presented mathematical analysis and show the tightness 
ofUne proposed approximations. 

, ., I. Introduction . _ 

Multihop communications, a viable option for providing 

broader and more efficient coverage, can be categorized 
as either non-regenerative (amplify-and-forward, AF) or re- 
generative decode-and-forward, DF) depending on the relay 
functionality |1|-|9|. In DF systems, each relay decodes its 
received signal and then re-transmits this decoded version. 
In AF systems, the relays just amplify and re-transmit their 
received signal. Furthermore, a system with AF relays can use 
channel state information (CSI)-assisted relays [l] or fixed- 
gain relays [2] (also known as blind or semi-blind relays Q). 
A (CSI)-assisted relay uses instantaneous CSI of the channel 
between the transmitting terminal and the receiving relay to 
adjust its gain, whereas a fixed-gain relay just amplifies its 
received signal by a fixed gain |2||6|. Systems with fixed-gain 
relays perform close to systems with (CSI)-assisted relays |2|, 
while their easy deployment and low complexity make them 
attractive from a practical point of view. 

Several works in the open literature have provided perfor- 
mance analysis of AF or DF systems in terms of bit error 
rate (BER) and outage probability under different assumptions 
of the amplifier gain ||T|-||9l. Among them, only two works 
dealt with the dynamic, time-varying nature of the underlying 
fading channel, |8|, |9|, despite the fact that it is necessary 
for the system's design or rigorous testing. In |[8l, the level 
crossing rate (LCR) and the average fade duration (AFD) of 
multihop DF communication systems over generalized fading 
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channels was studied, both for noise-limited and interference- 
limited systems, while Patel et. al in |9| provide useful exact 
analytical expressions for the AF channel's temporal statistical 
parameters such as the auto-correlation and the LCR. However, 
the approach presented in |91 is limited only to the dual-hop 
fixed-gain AF Rayleigh fading channel. 

In this paper, we study the second order statistics of the 
fixed-gain AF multihop Rayleigh fading channel. More specif- 
ically, motivated by the fact that this channel is a cascaded 
one, which can be modeled as the product of N fading 
amplitudes, we derive a novel analytical framework for the 
evaluation of the average LCR and the AFD of the product of 
N Rayleigh fading envelopes. Furthermore, we derive simple 
and efficient closed-form approximations using the multi- 
variate Laplace approximation theorem [16, Chapter IX. 5], 
iflTl . These important theoretical results are then applied to 
investigate the second order statistics of the multihop Rayleigh 
fading channel. Numerical and computer simulation examples 
verify the accuracy of presented mathematical analysis and 
show the tightness of the proposed approximations. 

II. Level Crossing Rate and Average Fade 
Duration of The Product of N Rayleigh Envelopes 

Let {Xi{t)}fLi be N independent and not necessarily 
identically distributed (i.n.i.d.) Rayleigh random processes, 
each distributed according to ifTOl - lfTTl . 
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in an arbitrary moment t, where ilj — E{Xf{t)} is the mean 
power of the i-th random process (1 < i < N). 

If {Xi{t)}fL-^ represent received signal envelopes in an 
isotropic scattering radio channel exposed to the Doppler 
Effect, they must be considered as time-correlated random 
processes with some resulting Doppler spectrum. This Doppler 
spectrum differs depending on whether fixed-to-mobile chan- 
nel ifTOl - llTTl or mobile-to-mobile channel llT2l - llT3l appears 
in the wireless communications system. In both cases, it was 
found that time derivative of i-th envelope is independent from 
the envelope itself, and follows the Gaussian PDF lfT0l - lfT3l 
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with variance calculated as 



Xi 
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If envelope Xi is formed on a fixed-to-mobile channel, 
then fi = /„jj where /„jj is the maximum Doppler frequency 
shift induced by the motion of the mobile station | lOj-LllJ. If 
envelope Xi is formed on a mobile-to-mobile channel, then 
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where and are the maximum Doppler frequency 
shifts induced by the motion of both mobile stations (i.e., the 
transmitting and the receiving stations, respectively) [fJI. It is 
important to underline that the maximum Doppler frequency 
in a fixed-to-mobile channel is /dmax = fmi, whereas the 
maximum Doppler frecjuency in a mobile-to-mobile channel 
is fd max = frm + /rm- The above results are essential in 
deriving the second-order statistical parameters of individual 
envelopes, as the LCR and the AFD |10l, |[TT1, CI- 

Below, we derive exact and approximate solutions for both 
of the above parameters for product of N Rayleigh envelopes. 
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We denote Y{t) as A^*Rayleigh random process or, at any 
given moment t, A^^Rayleigh random variable, following the 
definition given in llT4l . 



For some specified value {Xi} 



^ ^ {xi}fLi, the product 
X,. The LCR of y at 
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Y is fixed to the specific value y — Yii^i 
threshold y is defined as the rate at which the random process 
crosses level y in the negative direction |10|. To extract LCR, 
we need to determine the joint probability density function 
(PDF) between Y and Y, /yy-(y,2/)i and to apply the Rice's 
formula [11, Eq. (2.106)], 
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Our method does not require explicit determination of 
in order to determine analytically the LCR of the 
iV*Rayleigh random process, as presented below. 

First, we need to find the time derivative of (|5]l, which is 



Y = Y 
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Conditioning on the first — 1 envelopes {Xi}^^^^ = 
we have the conditional joint PDF Y and Y written 
as /y5>|^^...^^_^(2/,y|xi, ...,xjv-i). This conditional joint 
PDF can be averaged with respect to the joint PDF of the 
A^ — 1 envelopes {Xi}^Jl^ to produce the required joint PDF, 

/yy(y,2/) 



/ x-i—O J 

X fxt{xi) ■ ■ ■ fxM-i{^N~i)dxi ■ --(ixN-i (8) 

where to derive (O the mutual independence of the — 1 
envelopes is used. 



The conditional joint PDF /yy|_Yj...Xjv_i(2/: y^i, 
a;jv-i) can be further simplified by setting Y = y and using 
the total probability theorem, 

fYY\x,-x«^Sy^y\^^^ -^^N-i) 

= /Y|i'Xi-x„_i(yl2/'2;i,...,a;Ar^i) 

X /y|Xi---x„_i(2/|a;i, •■■,a;Ar-i) , (9) 

where each of the two multipliers in (|9]l can be determined 
from the above defined individual PDFs and their parameters. 

Based on (|7|i, the conditional PDF /y|yjfj...Xjv-i (^1?^' ^i' 
...jXjv-i) can be easily established to follow the Gaussian 
PDF with zero mean and variance 
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The conditional PDF of Y given {Xi\^S^^ = {xj,^^^ that 
appears in (|9]l is easily determined in terms of the PDF of the 
remaining A^-th envelope, 

/y|Xi---Xjv„i(y|a;i, ...,2:jv_i) 

/ N-1 \ N-1 

Introducing (|9]l and (fTTl i into (O, then ([8]l into (|6]l and changing 
the orders of the integration, we obtain 
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The bracketed integral in ( fT2] i is found using ( fTOb as 
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By substituting ([T]i and ( fTsT i into ( fT2] i. we obtain the exact 
formula for the LCR as 
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In principle, (fl4] i together with JTSl l provide an exact 
analytical expression for the LCR of the product of the product 
of TV Rayleigh envelopes (i.e., iV*Rayleigh random process 
lfT4l '). However, ( fT4T ) becomes computationally attractive only 
for small values of N, where it is possible to apply a numerical 
computation method (as Gaussian-Hermite quadrature). 

Note that, (fT4l l is transformed into a single integral when 
N — 2, which, after introducing ^ for i — 1,2 and changing 
integration variable x with new variable t according x — y/t, 
reduces to the known result [9, Eq. (17)]. 

The AFD of Y at threshold y is defined as the average time 
that the A^^Rayleigh random process remains below level y 
after crossing that level in the downward direction, 
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where Fy( ) denotes the cumulative distribution function 
(CDF) of Y . Fortunately, Fy{-) was derived recently in closed- 
form [14, Eq. (7)], as 

1 
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where G[-] is the Meijer's G-function [15, Eq. (9.301)]. 
A. An Approximate Solution for the LCR 

Next, we present a tight closed-form approximation of (fl4] i 
using the multivariate Laplace approximation theorem [16, 
Chapter IX. 5], ifTTl for the Laplace-type integral 

J(A) / m(x) exp(-A/i(x))dx, (18) 

where u and h are real-valued multivariate functions of x = 
[xi, - ■ ■ ,xn^i\, a is a real parameter and D is unbounded 
domain in the multidimensional space R^^^. 
A comparison of ( fT4l i and ( fTSl l yields 
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and A = 1. Note, that in the case of (fl4l l. all the applicability 
conditions of the theorem are fulfilled. Namely, within the 
domain of interest D, the function /i(x) has a single interior 
critical point x = [xi, • • • , xn-i], where 

,1/2 
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l/N ^ 



l<i< N~l, 
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which is obtained from solving the set of equations dh/dxi = 
0, where 1 < i < - 1. The Hessian (A^ - 1) x (A^ - 1) 
square matrix A, defined by [15, Eq. (14.314) ], is written as 

8/f)i ••• 4/V '»il7Ar_i 

8/1^2 ••• 4./^n2^N-l 
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By using induction, it is easy to determine that the A^ — 1 
eigenvalues of A are calculated as jii — 4/51^ for 1 < i < 
N ~ 2, and /ijv-i — 'IN/Q.^-i- Thus, all eigenvalues of A 
are positive, which, by definition, means that the matrix A 
is positive definite. By means of the second derivative test, 
since the Hessian matrix A is positive definite at point x, 
/i(x) attains a local minimum at this point (which in this case 
is the absolute minimum in the entire domain D). 
At this interior critical point x, 
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where ( l23T l is obtained using (O. Now, it is possible to 
approximate ( fTsT i for large A as 



(A'-l)/2 




X exp 



(25) 



It is well-know that the determinant of the square matrix is 
equal to the product of its eigenvalues. Hence, A can be 
written as 

A^22(^-i) r!ArAf22(^-i) 



det(A) 



rW-l 



(26) 



Although approximation ( |25] ) is proven for large A lfT6l - lfT7l , 
it is often applied when A is small and is observed to be very 
accurate as well. Similarly to [18|, we apply the theorem for 
A = 1. Therefore, the approximate closed-form solution for 
the LCR of A^*Rayleigh random process Y at threshold y is 
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(27) 



The numerical results presented in Section IV validate the 
high accuracy of the Laplace approximation applied for our 
particular case. 

Combining and (|27l) into (O, the AFD of the 

A^^Rayleigh random process Y at threshold y is approximated 
as 

V -1/2 
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III. Second Order Statiscs of Multihop 
Transmission 

Next, we apply the important theoretical result of the 
previous Section to analyze the second order statistics of the 
multihop relay fading channel. 
A. System Model 

We now consider a multihop wireless communications sys- 
tem, operating over i.n.i.d flat fading channels. Source station 
S communicates with destination station D through iV — 1 
relays Ti, T2,..., Tn~i, which act as intermediate stations from 
one hop to the next. These intermediate stations are employed 
with non-regenerative relays with fixed gain Gi given by 

with Go = 1 and Co = 1 for the source S. In ( |29] l, Wo^i is 
the variance of the Additive White Gaussian Noise (AWGN) 
at the output of the i-th relay, and d is a constant for the 
fixed gain d. 

Assume that terminal S is transmitting a signal s{t) with an 
average power normalized to unity. Then, the received signal 
at the first relay, Ti, at moment t, can be written as 

n{t)^ai{t)s{t)+wi{t), (30) 

where ai{t) is the fading amplitude between S and Ti, and 
wi{t) is the AWGN at the input of Ti with variance Wo,i. The 
signal ri is then multiplied by the gain Gi of the relay Ti and 
re-transmitted to relay T2 . Generally, the received signal at the 
fc-th relay Tk (k = 1, 2, ...,N - 1) is given by 



rfe(t) Gk-iak{t)rk-i{t) + Wk{t) 



(31) 



resulting in a total fading amplitude at the destination node 
D, given by 



N 



a{t) ^Yla^it)G^^l 



(32) 
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B. LCR and AFD of Multihop Transmissions 

If the fading amplitude received at node Ti, ai{t), is a 
time-correlated (due to mobility of Ti_i and/or Ti) Rayleigh 
random process, distributed according to ^ with mean power 
Cli — E{ai{t)}, then the i-th element of the product in 
(|32T i, Xi{t) = ai{t)Gi^i, is again a time-correlated Rayleigh 
random process, distributed according to ([T) with mean power 

Comparing (|5]l and (|32] |. we realize that the total fading 
amplitude at the destination station D (i.e., the received desired 
signal without the AWGN) is described as the A^*Rayleigh 
random process Y{t) = a{t), whose average LCR and AFD 
are determined in the previous Section. 

If all stations are assumed mobile with maximum Doppler 
frequency shifts fmS, fmO, /mi(l < « < for the source 

S, destination D and relays, respectively, then for the i-th hop 
fi = /m(4-i) + fnir with /^q = fmS and fmN = fmD, and 
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i=l 



f2 = f2 
J i J niS 
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Combining ( |27| i and ( [33] l, we obtain approximate solution for 
the average LCR of the total fading amplitude a at the output 
of a multihop non-regenerative relay transmission system. 



-, 1/2 




2/N 
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where $ is given by ( fTSl l. We see that (|34] | approximates the 
average LCR of the total fading amplitude for arbitrary mean 
power of the fading amplitudes f2j, arbitrary relay gains Gi 
and arbitrary maximal Doppler shifts 

Note that, for N — 2, (|34l i is an efficient closed-form 
alternative to the corresponding one [9, Eq. (17)] for the dual- 
hop case, which is shown in next section to be highly accurate. 

IV. Numerical Results and Discussion 

In this section, we provide some illustrative examples for 
the average LCR and AFD of the fading gain process of the 
received desired signal at the destination of the multihop non- 
regenerative relay transmission system model from Section III. 
The numeric examples obtained from the derived approximate 
solutions are validated by extensive Monte-Carlo simulations. 

We considered a multihop system consisted of a source 
terminal S, 4 relays, and a destination terminal D. The fixed- 
gain relays are assumed semi-blind with gains in Rayleigh 
fading channel calculated as [2, Eq. (15)] and [6, Eq. (19)] 



^2 1 /I 
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where 7^ = Cti/Wo,j is the mean SNR on the i-th hop, 
and r(-,-) is the incomplete Gamma function. Relay gain 
calculated according to ( l35b assures mean power consumption 
equal to that of a CSI-assisted relay, whose gain inverts the 
fading effect of the previous hop while limiting the output 
power at moments with deep fading. 

Depending on the stations' mobility, we used two different 
2D isotropic scattering models for the Rayleigh radio channel 
on each hop of the multihop transmission system. For the 
fixed-to-mobile channel (hop), we used the classic Jakes 
channel model ifTOll - lfTlll . For the mobile-to-mobile channel 
(hop), we used the Akki and Habber's channel model [12]- 
fT3|. The Monte-Carlo simulations of the latter were realized 
by using the sum-of-sinusoids method proposed in ||T9l - ll20l . 

More precisely, all mobile stations are assumed to induce 
same maximal Doppler shifts /,„, while the destination D is 
fixed. For all hops, (li = fl and Wqa — Wq- Thus, 7^ = 
7, Gi,sb — Gsb, and the mean of Rayleigh random process 
Xi{t) — ai{t)Gi-i,sb is calculated as 



U,i — exp 



r 0, 



= n ,2<i<N 



(36) 



whereas ili = il is selected independently from the AWGN, 
since Gq = 1. In this case, 

n N-1 



— h exp 



iV- 1 
7 



7 



(37) 
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Fig. 1. Average LCR, 7; = 7 = 5 dB 



Fig. 3. Average LCR, 7^ = 7 = 20 dB 
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malized threshold 



Fig. 2. AFD, 7, = 7 = 5 dB 



Fig. 4. AFD, 7, = 7 = 20 dB 



Note that, when introducing above scenario into ( l34b . a and 
Cl appear together as a/ V^. 

Figs. 1-4 depict the received signal's normalized LCR 
(Na/ fm) or normalized AFD (Tafm) versus the normalized 
threshold (a/V^) at 3 different stations along the multihop 
transmission system: at relay T2 (curve denoted by iV = 2), 
at relay T3 (curve denoted by iV = 3) and at the destination 
D (curve denoted by = 5). All comparative curves show 
an excellent match between the approximate solution and the 
Monte-Carlo simulations. 
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